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RECIPROCITY PRINCIPLE IN DUCT ACOUSTICS 

by Young-chung Cho 

• 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

Various reciprocity relations in duct acoustics have been derived on the 1 
of the spatial reciprocity principle implied in Green's functions for linear waves. 
The derivation includes the reciprocity relations between mode conversion coeffi- 
cients for reflection and transmission in nonuniform ducts, and the relation be- 
tween the radiation of a mode from an arbitrarily terminated duct and the absorp- 
tion of an externally incident plane wave by the duct. Such relations are well de- 
fined as long as the systems remain linear, regardless of acoustic properties of 
duct nonuniformities which cause the mode conversions. 

LIST OF SYMBOLS 
uniform duct elements (Fig. I) 

ambient sound speed 

radiation directivity function of j duct mode (far-field measure- 
ment at spherical angle (0,(p)) 

spherical Hankel function of the first kind 

I** 

spherical Bessel function 

|k| 

axial wave numbers of m and \l modes 
wave (or propagation) vector in free space 
k/k 


A, B 
c(x) 
fj $><P) 

(kr) 

i 

j £ <kr) 

k 

V k f. 

k 

A 

k 
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n 


unit normal vector outward from a surface 


q j,£m 

ft 

mn 


r 


constant matrix element 

mode conversion coefficient (first subscript is mode number of 
incident wave; second subscript is mode number of reflected 
wave) 

radial coordinate variable in a spherical coordinate system 


S 

S A» S B 

s £m «*,r) 


Sw 

^m/i 


surface of integration (Figs . 1 and 6) 

cross-sectional surfaces of duct elements A and B 
scattering matrix element for plane wave with incident angle (a , y) 
and spherical wave (£m) scattered by the duct 

surface made of duct wall 

mode conversion coefficient for transmission (first subscript is 
mode number of incident wave; second subscript is mode num- 
ber of transmitted wave) 


V 


volume 


x 


duct axial coordinate variable 


-x 


x o3 x' 
s’ s 

Y fm 


/5<V*s> 

V*s> 


spatial coordinate variables ((x^Xgi x^) in Cartesian coordinate 
system; (r, 8, in spherical coordinate system) 

spatial coordinate variables on a surface 

spherical harmonics (asterisk superscript denotes complex con- 
jugate) 

polar angle of -Tc 
complex function 

admittance of locally reacting boundary 



3 


7 

<5(x s ,%) 

9 , 6 ' 

P(x s ) 

V^y) 

# $ 
m n 

<P*<P' 


0) 

Subscripts: 


azimuthal angle of -k 

surface Dirac delta function 

polar angles in spherical coordinate system 

mean density of air 

amplitude of j duct mode induced by the external plane wave with 
incident angle (oi , y) 

normalized eigenfunctions 

azimuthal angles in spherical coordinate system 

acoustical velocity potential function (subscript A or B is used 
to indicate duct element if necessary) 

angular frequency of wave 


A, B 
a,b 


• uniform duct elements A, B 


v , r 


mode numbers (each represents, collectively, a pair of integers 
for a mode of a three-dimensional duct) 


Superscripts: 


a,b uniform duct elements A, B 

INTRODUCTION 

Modal analysis is a frequently used method of studying duct acoustics. The 
major merit of this method is that a mode remains undistorted in a uniform duct 
section. However, in the application of modal analysis to realistic systems, 
problems are often encountered. Some of the difficulties are due to duct nonuni- 
formities such as duct termination, variation of duct cross section and wall im- 
pedance, alteration of duct axis (duct bending), 'axial gradients of mean fluid var- 
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iables, etc. The complexities can be reduced by virtue of basic physics laws 
such as the reciprocity principle. Furthermore, such a principle helps one to 
gain physical insights into the problems, often without detailed analyses. This 
paper discusses in some detail the reciprocity principle in duet acoustics. 

The reciprocity principle implies a symmetry property between a source and 
the excited field in systems of linear vibration . 1 The symmetry property is well 
elucidated in Green's functions for waves. The reciprocity principle applies to 
a variety of systems of linear vibration, including those which are bounded, un- 
bounded, or partially bounded. A ducted acoustics system may be considered as 
one which is partially bounded. Because a physics law should not depend on rep- 
resentations which are chosen for description of the system, the reciprocity prin- 

10 

ciple, which is well established in spatial coordinates, should be equally well 
deseribable in eigenfunction spaces. Note that eigenfunctions in a duct are in one- 
to-one correspondence with the duct modes. 

Some limited discussions of the reciprocity principle with the duct mode rep- 
resentation are given in a number of previously reported works. Levine and 
Schwinger briefly discussed the reciprocity relation between the radiation direc- 
tivity of the fundamental mode and the absorption of a plane wave by the duet.^ 

The reciprocity relation between mode conversion coefficients for reflection of 
the propagating duct modes was used by Weinstein for discussion of the Wiener- 

Hopf solutions of radiation problems from hard-walled uniform circular and 

5 

two-dimensional wave guides. In his numerical studies of sound transmission 
through two-dimensional duct bends, Tam found that the mode conversion coeffi- 
cients satisfy the reciprocity relation. Eversman discussed the reciprocity 
principle for sound reflection and transmission in a two-dimensional hard-walled 
duct with partially variable cross section . 1 Koch found that the reciprocity rela- 
tion is also satisfied by numerical results of the mode conversion coefficients for 
reflection and transmission in partially lined rectangular ducts of uniform cross 
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O 

section. A more fundamental view of the reciprocity principle in duct acoustics 

n 

was offered by Mohring. His approach is not, however, applicable to a lossy 
system like a lined duct, because the scattering matrix fails to be unitary in such 
a system. 

As noted in Ref. 1, Chapter V, the reciprocity principle applies to linear vi- 
bration systems involving dissipative force fields as well as conserved force 
fields if the dissipative forces are linear functions of the particle velocity. A 
wide range of dissipative forces belong to this type, including those pertinent to 
lossy boundaries involved in acoustic liners of ducts. 10 The reciprocity principle 

Q 

also applies to systems involving an inhomogeneous medium, inclusive of gra- 
dients of mean density, pressure, and temperature, as long as the systems re- 
main linear. 

In this paper, the reciprocity relations are derived for mode conversions 
caused by the duct nonuniformities which have been mentioned thus far. The de- 
rivation is based on the reciprocity principle implied in Green’s functions for 
waves, and includes the reciprocity relation between mode conversion coefficients 
for reflection and transmission in nonuniform ducts and for reflection due to an 
arbitrary duct termination. Also included is the reciprocity relation between the 
radiation directivity of a duct mode and the excitation of the mode from the ab- 
sorption of an externally incident plane wave by the duct. 

RECIPROCITY RELATIONS IN NONUNIFOEM DUCTS 

In this section, we consider sound reflection and transmission in a nonuniform 
duct as schematically illustrated in Fig. 1. The duct is composed of two uniform 
duct elements A and B, which are connected through a nonuniform transition sec- 
tion. The two uniform elements may differ from each other in cross-sectional 
shape and area, duct axis, wall impedance, and mean fluid variables. Also shown 
in the figure are the surfaces s A , s B , and S^., which are, respectively, the 
cross sections of duct elements A and B, and the duct wall. 
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On the assumption of harmonic time dependence e - *^, the wave equation is, 
in the absence of mean flow, reduced to 

v 2 m> (x) + (x) = 0 (1) 

where the sound speed c(x) may vary due to a spatial gradient of the mean fluid 
variables . The analysis includes boundary conditions on the duct walls which can 
be given by 

n • VT(x s ) + / j3(x s ,xp* (xpd 2 x^ = 0 (2) 

* /S W 

where /3 is a complex function of x_ and xL. This relation is a general form of 

S S 

homogeneous boundary condition. Note that boundary conditions for eigenvalue 

9 

problems are always homogeneous. This equation can be used to account for 
duct walls of extended reaction as well as for locally reacting walls. For the 
latter, /3 is given by 

/3(x s ,3?) = -iajp(xp/? o (xp6(x s ,x^) (3) 

where /3 q is the wall admittance, and 6(x s ,x^) is the Dirac delta junction. 

Reciprocity for Transmission 

We first derive the reciprocity relation for sound transmission in the duct. 

We use a method similar to one which is used for the proof of the spatial reci- 
procity principle in Ref. 3. Consider two different incident waves. One is com- 
posed of a single mode m generated In the duct element A (see Fig. 2), and the 
other of a single mode fi generated in the duct element B (see Fig. 3) . Both in- 
cident waves propagate towards the nonuniform section. 

The incident wave from A, upon arrival at the nonuniform section, is partly 
reflected back to A and partly transmitted to B, as schematically illustrated in 
Fig. 2. The wave will be attenuated unless the duct wall is hard. The resultant 
wave, denoted by ^ A , can be written in the uniform duct elements as follows: 
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Transmitted from 


B to A 


Note that mode numbers m, n, and are used for duct element A; and /i, v, 
and r are used for element B. Each mode number represents, collectively, a 
pair of integers 'for a mode of a three-dimensional duct. The sub- or superscript 
a or b is used on $ , k, x, and R to indicate the corresponding duct element, 
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A or B. For the transmission coefficient T, the superscript refers to the duct 
element in which the incident wave is generated. It should also be noted that the 
two sets of eigenfunctions, and can differ from each other. 


Let us consider the volume integral 


I = / 


( 6 ) 


where V is the volume enclosed by the surfaces S A’ S B’ and S w . It follows 
from Eq. (1) that the integrand is identically zero everywhere. Thus, the inte- 
gral is equal to zero. By virtue of Green's theorem, the integration is trans- 
formed into 



( 7 ) 


It is convenient to carry out the integration separately on the surfaces S^, Sg, 
and S^y. It can be readily shown by using Eq. (2) that the integration over S w 
vanishes. The integration on is written as 



( 8 ) 


Substituting Eqs. (4) and (5) into Eq. (8) and utilizing the orthonormality of eigen- 
functions yields 


*A - 

Similarly, the integration on Sg results in 

Ir> = 2iT?_ 

-b ttl/I fl 

It follows from Eqs. (9), (10), and (7) that 

T 3 - Ir^ _ 'T’fr 
m fi ij, fim m 


(9) 

(L0> 


( 11 ) 
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This equation specifies the reciprocity relation between the mode conversion coef- 
ficients for transmission, or simply, the transmission coefficients. Recall that 
the transmission coefficient is the amplitude of mode p of duct element B, 

which is induced by the incidence of mode m of duct element A. On the other 
hand, T^ m is the amplitude of mode m of duct element A , induced by the inci- 
dence of mode p of duct element B. Once the value of is computed, then 

the value of T* 3 is determined (or vice versa) from Eq. (11) . 

ptm 

Reciprocity for Reflection 

For the reciprocity relation for reflection, we consider two waves, ^ and 
Wave is the resultant wave due to the incidence of a new single mode 
j (^m) of duct element A. It can be written in the uniform duct elements as 
follows: 

At A, 


At B, 


a 

a ik i x a 
= 3 a + 


Incident 


Z a a ~ ik n x ' 

E ?n*n e nc 


-L 


v 

Reflected 


iiJb. 


n- ? 


i t ; 

Transmitted from 
A to B 

where j ^ m. We then have, in the place of Eq. (7), 


(12a) 


(12b) 



( 13 ) 
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The integral over again vanishes. With the substitution of Eqs. (4) and 
(12) into Eq. (13), the integration on S ^ and Sg is easily carried out. The con- 
tribution from the integration over Sg is also zero, and the integration over 
gives 

CL4) 


It follows from Eqs. (13) and (14) that 

B®,k? = Rf_k* (15) 

mj j jm m 

This is the reciprocity relation for reflection. Note that this reciprocity relation 
is valid for reflection due to any type of duct nonuniformity if the linearity re- 
quirement is fulfilled, that is, if Eqs. (1) and (2) are valid, respectively, in the 
media and on the boundaries . This can be readily proven from the fact that in 
the derivation of Eq. (15), the integrations on Sg and make no contributions 
to the integral in Eq. (13). 

It should be mentioned that the reciprocity relations derived here are valid 
for the whole frequency range. In other words, Eqs. (11) and (15) are valid for 
real, imaginary, or complex values of axial wave number. 

RECIPROCITY RELATION IN RADIATION PROBLEM 


In this section, we derive the reciprocity relation between the radiation di- 
rectivity of a duct mode and the excitation of the mode due to the absorption of an 
externally incident plane wave. The geometry of the problem is schematically 
illustrated in Figs . 4 and 5 . 

We consider again two incident waves. One is composed of the j mode of the 
uniform duct element A, which is propagating toward the duct termination. This 
wave is partly reflected back to duct element A, and is partly radiated. It will 
be also partly attenuated if the duct wall is soft. The resultant wave, denoted by 
can be expressed as follows (see Fig. 4): 
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In duct element A, 


ik-x V - ' -ik x 

*A =4 i e 


n 


Incident Reflected in duct 
in duct 


(16a) 


where the super- or subscript a has been dropped from the variables without the 
loss of clarity. The radiated wave can be written in terms of the spherical Hankel 
function of the first kind and the spherical harmonic, the product of which will be 
referred to as outgoing spherical wave functions : 

’ t A = E q j.<!m h i 1> < kr ) Y fm( 9 >' ! ’> < 16b > 

Vim ; I / 

Radiated from duct 

Here, the radiation is into the three-dimensional space; l and m are single in- 
teger numbers; and qj is a constant matrix element. 

In the far-field region, Eq. (16b) can be written as 


where 




f.(0,<p) 


(17) 


¥ 8 ’ f ’> = # ( 18 ) 

fm 

It follows from Eq. (17) that f- (0, cp) completely determines the radiation direc- 
tivity of the duct mode. It will be herein referred to as the radiation directivity 
function. 

The other incident wave under consideration is an externally incident plane 
ik * x 

wave e . Since the present problem deals with a constant sound frequency, 
the magnitude of the wave (or propagation) vector is predetermined by k = (cu /c) . 
Thus, this incident wave is specified primarily by the direction of its propagation. 
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k. We will use the incident spherical angle (a, y), which is in one-to-one corre- 

A 

spondence with k. 

The externally incident wave is partly scattered by the duct and partly trans- 
mitted into the duct, as shown in Fig. 5. The scattered wave can be expressed in 
terms of the outgoing spherical wave functions . The resultant wave outside the 
duct is then written as 


*B = 7=r eik ' X+ Z 

J 2ir / j 


(19a) 


Incident Scattered 

s 

Here, the constant of the first term is the normalization factor of the plane wave 
in the sense of the Dirac delta, and £>£ m ( a > y) is the scattering matrix element. 
The wave transmitted into the duct element A ^Is written as 


* 


B 




(a,Y)$ n e 


-ik x 
n 


, 1 *i 


Transmitted into 
.duct A 


(19b) 


where o n (a , y) is a constant matrix element. 
We now consider the integral 




*A W B~*B W A 


)•* 


dS = 0 


( 20 ) 


where the surface S is shown in Fig. 6. This integral vanishes as long as 
Eq. (1) is valid throughout the medium enclosed by S. The contribution from the 
integration on the inner and outer walls of the duct is zero if the wall boundary 
condition can be given as in Eq. (2), Consequently, we have 


( 21 ) 
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where 



( 22 ) 



(23) 


Inserting Eqs. (16) and (19) into Eqs. (22) and (23) and utilizing Eqs . (18) and (21) 
yields 

27rfj (a,y) = ik. oj (a , y) (24) 

Here, in evaluating the integral I B , we have used the various function properties 
given in Appendix A. 

Equation (24) specifies the reciprocal relation between the radiation direc- 
tivity of a duct mode and the amplitude of that mode induced by external plane 
waves. Note that oj (a , y) is the amplitude of the j mode of the duct, which is 
induced by an external plane wave incident to the duct termination with the incident 
angle (o-,y). On the other hand, f - (a , y) is the amplitude of the far-field measured 
at the spherical angle (a ,y) for the radiation of the j mode. If the radiation di- 
rectivity of a duct mode is known, the amplitude of the mode induced by an exter- 
nal plane wave is computed by this equation. Or conversely, the radiation direc- 
tivity of a duct mode can be indirectly determined by measuring the amplitude of 
the mode induced by external plane waves of various incident angles. 

CONCLUDING REMARKS 

Various reciprocity relations with modal analysis of duct acoustics have been 
derived. The existence of such relations are solely due to the linear characteris- 
tic of the vibration systems under consideration. The derivation is general and 
is valid regardless of the complexities of the duct nonuniformities involved. The 
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reciprocity relations are-useful in studies of sound propagation in nonuniform 
ducts and radiation from a duct. For instance, the numerical computation in- 
volved in such studies can be reduced by means of the reciprocity relation. For a 
radiation problem, the far-fieid directivity can be studied indirectly by measuring 
the transmission of external plane waves into the duct. Finally, it is always ad- 
vantageous for studies of a system to have a rule which is not affected by compli- 
cations of the system. 
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APPENDIX A 

EQUATIONS PERTAINING TO SPHERICAL WAVES 
Listed here are the various equations used in deriving Eq. (24). 

The orthonormality of the spherical harmonics: 

sin 6 de d( P ~ 6 W 6 mm> < Al ) 

The completeness of the spherical harmonics: 

co 1 

V V ? &m (d,<p)Yl m (e\<p') = 6{cose-.cose')5(<p-<p') (A2) 

£=0 m=-£ 

The relation between the plane wave and the spherical harmonics: 

«> i 

e ik.x _ 4 ff ^ (A3) 

£.= 0 m=-£ 

where the vectors are given in terms of spherical coordinate variables as (see 
Fig. 5) 

x = (r, 6, <p) 
k = (k,a! , ,y) 

and jjj (kr) is the spherical Bessel function. 

Also used is 

Y i m ( d > ( P) = (~ 1 ) &Y n m (' n ~ 7r+ <P) ( A4 ) 

Note that the direction corresponding to the spherical angle (ir - 0, it + <p) is the 
reverse of the direction corresponding to the spherical angle (0,<p). 
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The asymptotic form of the spherical Hankel function: 



h / +1 e 


ikr 


•kr 


(A5) 
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Figure 1. - Schematic of non uniform duct. 
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Figure 2. - Schematic representation of wave function 
% given in Eqs. (4a).and (4b). 
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Figure 3. - Schematic representation of wave function. 
'J'b given in Eqs. (5a) and (5b). 



Figure 4. - Sound radiation from an arbitrarily terminated duct. 
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Figure 5. - Scattering and absorption of an external 
plane wave by a duct. 



Figure 6. - Surface of integration, S, in Eq. <20). 
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Figure 7. - Example case of sound reflection due to 
' duct nonuniformity. 
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